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ABSTRACT 

We conducted micro-gravity experiments to study the outcome of collisions between sub-mm-sized 
dust agglomerates consisting of /im-sized Si02 monomer grains at velocities of several cms~^. Prior 
to the experiments, we used X-ray computer tomography (nano-CT) imaging to study the internal 
structure of these dust agglomerates and found no rim compaction so that their collision behavior is 
not governed by preparation-caused artefacts. We found that collisions between these dust aggregates 
can lead either to sticking or to bouncing, depending mostly on the impact velocity. While previ- 
ous collision models derived the transition between both regimes from contact physics, we used the 
available empirical data from these and earlier experiments to derive a power law relation between 
dust-aggregate mass and impact velocity for the threshold between the two collision outcomes. In 
agreement with earlier experiments, we show that the transition betwe en both refflmes is not sharp, 
but follows a shallower power law than predicted by previous models (jGiittler et al.|[2QTQ ). Further- 



more, we find that sticking between dust aggregates can lead to the formation of larger structures. 
Collisions between aggregates-of- aggregates can lead to growth at higher velocities than homogeneous 
dust agglomerates. 



1. INTRODUCTION 

Within the past years, close collaborations between ex- 
perimental and theoretical physicists enriched our under- 
standing of the processes which finally lead to the for- 
mation of planet esimals. From the experimenters' side 
many experiments o n the collisions behavio r (se e the ex- 



perinient rev iews by Blum & Wurm ( 2QQ8t and Guttler 
et al. (2010)) and mechanical properties ( [Guttler et al. 
2009p of porous dust agglomerates were conducted. Lab- 
oratory experiments have shown that collisions can lead 
to various outcomes. Dust agglomerates will stick, due 
to the inelasticity of the collisions and the attractive 
van-der-Waals force, as long as their collision velocity 
does not exceed the bouncing threshold velocity. Stick- 
ing can occur either at the first contact at very low ve- 
locities or after some plastic deformation in the contact 
region. Guttler et al.| (2010) developed an experiment- 
based dust- aggregate collision model, which predicts the 
outcome of collisions between dust agglomerates as a 
function of the masses, porosities, and collision veloc- 
ity of t he du st aggregates. Based on this model, Zsom 
et al. (2010) simulated the growth of dust aggregates 
m various protoplanetary-disk (PPD) models, usin^ a 
Monte-Carlo code. The simulations by Zsom et al. (2010) 
showed that, depending on the PPD model, ratEer com- 
pact dust aggregates of sizes up to one centimeter can 
grow at 1 AU. However, the growth always stops due to 
the occurrence of bouncing in dust- aggregate collisions, 
which reduces the porosity and does not lead to a further 
increase in aggregate mass. This has been re ferred to as 
the bounc ing barrier. More recent work by Windmark 
et al. ( 2012a b ) showed that the bouncing barrier can be 
overcome either by the introduction of a few bigger par- 



ticles or by taking a Maxwellian distribution of collision 
velocities into account. Both processes lead to the onset 
of a fragmentation-induced mass transfer from smaller to 
larger bodies. 

However, it has to be noted that there were only a 
few experiments which studied collisions between dust 
aggregates close to or at the threshold between sticking 
and bouncing so that all previous outcomes are somehow 
model depende nt and not empirically proven. Recently, 
Weidling et al.' (2012) (in the following referred to as Pa- 
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per I) studied collisions between mm-sized dust agglom 
erates. They found that the transition between sticking 
and bouncing was not sharp, but followed a wide transi- 
tion from perfect sticking to perfect bouncing, spanning 
a range of two orders of magnitude in velocity. As the 
impact velocities in the experiments by Weidling e t al.| 
(2012) were not sufficiently low to study the whole tran- 
sition range, they modified a collision mo del for plas- 
tic/elastic spheres ( Thornton fc Ning||1998 ) for dust ag- 
gregates, taking into account the porous surfaces of the 
colliding bodies. They confirmed the original thresh- 
old between sticking and bouncing introduced by Gi ittler| 
|et al.| ( |2010| ) as the line where half of the dust aggregates 
stick. 

To gain a deeper insight into the physics of dust- 
aggregate collisions, numerical simulations have made 
increasing contributions. Internal processes, like restruc- 
turing and the formation of cracks in a collision, are not 
easily observable in laboratory experiments but are read- 
ily available in numerical simulations. They also allow 
for the study of collisions of dust agglomerates with in- 
ternal structures (e.g., very high porosities, well-defined 
inhomogenities) that are hardly feasible in experiments. 
However, one major discordance between experiments 
and simulations has not been clarified so far: in most 
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cases simulations are not able to achieve bouncing in the 
way it was observed in several experime nts (HeiBelmann 



et al.||2QQ7 Weidling et al.||2Q12 



Blum fc Miinch|,1993J . 



in these cases, collisions between dust aggregates lead 
to the rebound of the particles, without any visible mass 
lo ss, mass transfer, or d e forma tion of the dust aggregates. 



Geretshauser et al. 



^ ^ ( 2Qll| compare d the collision 

outcom es observed and categorized by [Guttler et al. 
(2010) with the results from their smoothed parti- 
cle hydrodynamics (SPH) model, w hich had been cali- 
bra ted w ith laboratory experiments dGiittler et al.||2009 



Geretshau ser et al.||2010| ). [Geretshauser et al.~i 2nTT 
were able to observe bouncing collisions in their simu- 
lations. They collided spherical dust agglomerates with 
10 cm and 6 cm diameter and an initial volume filling 
factor of ^ = 0.35 (i.e., with a porosity of 65%) at differ- 
ent velocities. The simulations showed that the particles 
bounce at velocities of less than 1 ms~^. The colliding 
dust aggregates are visibly compacted to ^ = 0.45 within 
the contact volume. At higher impact velocities, suffi- 
cient energy is dissipated in restructuring so that the con- 
tact area between the two spheres becomes large enough 
to lead to sticking. This is in contrast to th e experimen- 



tal observations collected in [Guttler et al. (2010) where 
sticking only occurs at (mucn) lower velocities. Sticking 
in collisions between dust aggregates at high velocities 
was experimentally only observed if one of the dust ag- 
gregates was much smaller (projectile) than the other one 
(target). In these cases, the deep penetration of the pro- 
jectile into the ta rget agglomerate leads to the capturing 
of the projectile ( Langkowski et aT[|2008 ). Experiments 
with cm-sized or larger dust agglomerates with various 
porosities have shown no sign of stic king at or ar ound 
1 ms-^ ("Beitz e t alT1|2Qri[ [Schrapler et al.[[2012[ ). In 
contrast, all laboratory experiments have contirmed the 
transition between bouncing and fragment ation at a ve- 



locity of about 1ms veloc ity (.G iittler et al. 
et al.|[2011l [Schrapler et'll|2012| r 
Ke bound in central collisions ha 
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Beitz 



"has also been stu died in 
the numerical simulations by Wada et al.[ (|2011|). They 
used a soft-sphere molecular-dynamics model to simu- 
late collisions between dust agglomerates consisting of 
^ 4000 monomer grains with various filling factors and 
collision velocities. The monomers were spherical Si02 
particles with a diameter of 0.1 /im. The porosity of 
dust agglomerates can be described by the coordination 
number, which is the average number of (touching) next 
neighbors of a monomer grain. For the case of a close- 
packed hexagonal structure, from w hich a fraction fex 



of particles are randomly withdrawn, Wada et al. ( [2011 ) 
give the following dependence of the volume tilling factor 
on the coordination number Uci 
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where (j) is the volume filling factor of the dust aggregate 
and (f)cp = 0.74 is the filling factor of the hexagonal-close 
packing. They also studied the rebound behavior of cubic 
latice structures and agglomerates produced by ballistic 
agglomeratio n with one (B AMl) and two (BAM2) mi- 
grations (She n et al.[|2008 ) made out of ice particles. It 
is difficult to compare the former, highly arranged cu- 
bic particle with natural grown agglomerates. The latter 



agglomerates are of interest. The BAMl and BAM2 ag- 
glomerates have coordination numbers of 4 and 6 and 
volume filling factors of ^ 0.23 and ~ 0.34, respectively. 
The simulatio ns for agg lomerates, based on hexagonal 
packing, by Wada et al. ( [20^1 1| ) have shown that highly 
porous dust aggregates (low coordination number) stick 
to each other even if their collision velocity is 22 ms~^ 
(mind that the discrepancy between this high value and 
the fragmentation threshold of '^1 ms~^ is mostly due 
to the factor of 15 in monomer diameter and the usage 
of two different materials, i.e. Si02 and H2O ice). Only 
if the coordination number reaches values larger than 
6, the dust aggreg ates (based on closest packing) in 



the simulations by Wada et al. (2011 ) bounce off at inter 
mediate (2.2 ms~^) and high velocities (22 ms~^). With 
Eq. [ij this coordination number corresponds to a volume 
filling factor of approximately (j) = 0.37. In the studies of 
the two BAM agglomerates, bouncing was also found for 
BAM2 agglomerates with nc = 6, only. For these aggre- 
gates no dependency on the velocity was found. |Wada] 
et al. (2011) also found that this critical value for the 
coordination number is independent of the agglomerate 
material, because above these values the particles are me- 
chanically locked. This prohibits the dissipation of the 
collision energy by r estructuring of the agglomerate and 
leads to bouncing. |Wada et al.| (|2011|) argue that the 
bouncing observed in experiments might be caused by a 
preparation-induced compacted surface layer or the over- 



all fill ing factor of the agglomerates. While |Wada et aL 
(2011) can measure the volume filling factor and the co 



ordination number of their simulated aggregates, we can 
only measure the filling factor and, therefore, depend on 
a relation between these two values. The relation fo r the 
coordination number given by Wad a et al.[ ( [201l| (Eq. 
[1]) yields higher coordination numbers than relations for 
oifferent structures referred to in other publications. An 
overview a bout some other relations is given in Fig. [l] 
Van de Lagemaat et al. (2001) calculated the dependence 



of the coordination number on the volume filling factor 
for a randomly packed Ti02 nano-particle film and found 
a coordination number of tt-c = 6 for a volume filling fac- 
tor of (/) = 0.57. Furthermore, [Van Antwerpen et al. 



(2010) gave a summary of different simulations and ex- 



periments and according to this paper, the highest coor- 
dination numbers (in the range releva nt to this work ) are 
achieved by the equation deriv ed by [Rumpf ( 1958| ) and 
the lowest by 'Yang et al. (2000). Following this, the crit- 
ical volume filling factors lor bouncing (i.e., correspond- 
ing to Tie = 6) are (j) = 0.48 and (j) = 0.65, respectively. 
Recent numerical studies by Schmidt and Blum (pers. 
comm.), who simulated the growth of agglomerates by 
random ballistic deposition (RBD), where the particles 
stick at the first contact, and random gravitational de- 
position (RGD), where particles were allowed to roll and 
bounce, ha ve also yielded compara ble results (triangles 
in Fig. [1]). [Seizinger fc Kley (2013) simulated collisions 
between microscopic dust aggregates. To produce ag- 
glomerates which are as realistic as possible they also 
used a ballistic aggregation with migration method. An 
agglomerate is produced by monomers impacting from 
random directions on the target. After the first contact 
with the agglomerate the monomer can either stick or 
migrate to a position where it is in contact with two 
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Fig. 1. — Comparison of different numerical models describing 
the relation between coordination number and volume filling fac- 
tor of packed particles. The triangles represent the simulation from 
Schmidt and Blum (pers. comm.) for random ballistic deposi- 
tion (RBD) and random gravitational deposition (RGD). The dot- 
ted vertical line represents the average volume filling factor of our 
dust agglomerates and the shaded range denotes the ±lcr range. 
The squares represent the coordination number and filling factor of 
different lattices (diamond, simple (SC), body-centered (bcc) and 
face-centered (fee) cubic, resp. hexagonal closed packing (hep)). 



or three monomers. The selection of the final position 
has been chosen by three different ways: By selecting 
the closest position to the impact point (BAM short- 
est), by selecting a random position (BAM random), or 
by selecting the position which is closest to the center 
of the agglomerate (BAM center). As shown in Fig. [ij 
the three methods lead to different coordination numbers 
for different volume filling factors. All of these methods 
lead to coordination numbers, comparable to the other 
models, but l ower t han the BAM agglomerates used by 
Wada et al. (2011). This is caused by different calcu- 
lation methods of the filling factor. Seizinger fc Kley| 
([2013') used BAM agglomerates which were 
of a larger agglomerate. 



sliced out" 
Therefore, irregularities on the 
agglomerate's surface do not decrease the filling factor. 
Also their agglomerat e s are considerable bigger than the 
ones from Shen et al. (2008), which also reduces surface 
effects. This makes the agglomerat es more comparable to 
the particles used, for example, in'Weidling et al.|(|2012|) 
or this paper. In this paper, we will show that bouncing 
in aggregate- aggregate collisions occurs for lower coor- 
dinati on numbers than those predicted by |Wada et al. 
(I20TTI). 



In Section |2.1 1 we will present our experimental setup, 
which allows us to study low-velocity collisions between 
dust agglomerates and their growth under micro-gravity 
conditions (also see Paper I). We will discuss the prop- 
erties of the experimental d ust a gglomerates in Section 
|2.2[ Subsequently, in Section [231 we will describe the in- 
ternal structure of our dust aggregates with respect to 
the discussion above. The results of our collision experi- 
ments will be described in Section [3] where we will focus 
on the transition between sticking and bouncing and the 
formation of larger dust clusters. Section |4l deals with 
the modeling of the data, and a conclusion of the paper 
is given in Section [5j 



2. EXPERIMENTAL SETUP 

In this Section, we will briefiy describe the experimen- 
tal setup and the dust agglomerates used in our study. 
A more detailed description of the experimental method 
can be found in Paper I. 

2.1. Multiple Particle Experiment 

Studies of collisions between dust aggregates at very 
low impact velocities are difficult to achieve in the labo- 
ratory due to the presence of gravity. Therefore, we used 
the Bremen drop tower to generate a large number of 
collisions between dust agglomerates under micro-gravity 
and vacuum conditions. The setup of this experiment is 
based on the shaking device originally introduced in Pa- 
per I (see Fig. [2]F). It consists of a particle container, a 
shaking mechanism to agitate it and thus generate the 
collisions, an LED array with diffusor screens for back- 
light illumination, and a high-speed camera to record 
the dust- aggregate motion and collisions during the ex- 
perimental runs. In addition to the setup described in 
Paper I, the much smaller particles are contained in a 
small cylindrical glass housing with 1 cm diameter and 
a height of 1.6 cm. The glass cylinder itself is mounted 
inside a larger glass vacuum chamber, which is connected 
to and thus evacuated by the drop tower . To reduce the 
pressure even more we additionally used a turbomolecu- 
lar pump. A fine-meshed sieve was attached to the top of 
the inset container to act as a wall for the particles but 
enable us to evacuate the smaller volume. As described 
in Paper I, the vacuum chamber is attached to a platform 
which can be shaken by an eccentric wheel driven by a 
motor. With a fixed amplitude of the wheel, the velocity 
with which the particles are agitated inside the vacuum 
chamber is changed with the rotation frequency of the 
motor. Once the experiment is under micro-gravity con- 
ditions, the shaker distributes the particles through the 
container (see Fig. [23). The dust aggregates are illu- 
minated by a back-light LED array so that the dust ag- 
gregates appear dark in front of a brighter background. 
To record the experiment, we used a high-speed camera 
with a resolution of 512 x 512 pixels and a repetition rate 
of 500 frames per second. The camera observes the ex- 
periment through a prism, which generates two images 
separated by a viewing angle of 30 degrees. 

In this paper, we analyzed four experimental runs of 
the experiment (see Table [T]). We used two dust-analog 
materials described in the next section in two experi- 
ments each. The total duration of each of the four ex- 
perimental runs was ~9 seconds. For experiments 1 and 
2, the shaking frequency was reduced to 50% of the ini- 
tial value after 5 seconds. In experiments 3 and 4, the 
frequency was reduced to 50% of the initial value after 
2 seconds, to 35% after 5 seconds, and the motor was 
completely shut off after 8 seconds. 

2.2. Characterization of the Dust Aggregates 

The experiments were conducted with two different 
dust analog materials as shown in Fig. |2^ and [2)3. The 
first analog material is the same as used in Paper ^ It 
consists of Si02 particles with an irregular shape and a 
wide size distribution, shown in Fig. [3j The average 

^ Manufacturer: Sigma- Aldrich 
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Fig. 2. — a. SEM image of polydisperse, irregular dust grains with a size distribution of 0.5 - 10 /im. b. SEM image of monodisperse, 
spherical dust particles with a diameter of 1.5 /im. c. Computer tomography image of a dust agglomerate consisting of monodisperse Si02 
particles prior to the experiment, d. Optical microscopy image of sieved dust agglomerates, e. Snapshot from one of the experiments, f. 
Experimental setup (adapted from Paper I). 



TABLE 1 

Overview of the performed experiments. The colors in the time table indicate the PERroRMANC^)F the shaking mechanism 

FROM THE START OF THE EXPERIMENT (SoE) TO THE END OF THE EXPERIMENT (EoE). ^^^^ , , 35% , AND 0% . 
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particle size by number is 0.63 /im and by mass 2.05 
jam. Approximately 80 % of the particles are smaller 
than 1 /im, while more than 80 % of the mass is in par- 
ticles in excess of 1 /am diameter. The mass density of 
the monomer grains is 2.6 gcm~^. The second material 
consists of Si02 spheres with a monodisperse size distri- 
bution and a diameter of 1.5 /inr] The mass density of 
this material is 2.0 gcm~^. Both materials have been 

^ Manufacturer: micromod Part ike Itechnologie GmbH 



used in many earli er experiments and are listed as 
and Qf4 in Blum & Wurm (2008). 

The dust agglomerates used in our experiments were 
formed in the dust-storage containers. To prepare the 
dust agglomerates, we sifted the dust through two sieves 
with a mesh size of 250 /im (upper sieve) and 100 /im 
(lower sieve). The selected dust agglomerates have an 
average size of 150 ±60 /im(standard deviation) and have 
a roughly spherical or ellipsoidal shape (see Fig. |2]i). 

The volume filling factor of the sieved dust aggregates 
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( 2000 1) , the measured volume filling factor rather corre- 



Particle Size [pm] 

Fig. 3. — Cumulative size distribution of the irregularly shaped 
Si02 particles. Red squares: cumulative frequency by count of 
the particles. Blue spheres: cumulative frequency by mass of the 
particles. 

was derived in Paper I for several specimen by measure- 
ments of the dust- aggregate masses with a high preci- 
sion balance and their respective volumes with optical 
methods. In Paper I, an average volume filling factor 
of (/) = 0.35 ± 0.05 was determined for mm-sized dust 
aggregates with masses between 0.1 and 30 mg. The 
dust aggregates studied in this paper are considerably 
smaller so that the determination of the mass and vol- 
ume of individual dust aggregates was not feasible with 
this method. 

However, to get a more accurate measurement of the 
volume filling factor and also to study the inner struc- 
ture of our dust aggregates (for those with monodisperse 
monomers), we used nano-CT images obtained with a 
SkyScan 2011 scanner with a voxel size of 0.3 (/im)^. We 
analyzed one dust aggregate with a triaxial ellipsoidal 
shape and semi-major axes of approx. 150 /im, 100 /im, 
and 40/im, prior to the experiments. We calculated the 
volume filling factor by measuring the average gray value 
in the images and normalized this to the maximum value 
in the sample, which represents the density of solid Si02. 
The difficulty of this method is to distinguish between the 
gray-value of the bulk material and artifacts caused by 
the reconstruction of the tomography images. To calcu- 
late the gray value for Si02, we measured the maximum 
gray values inside of 50 monomer grains. At the given 
resolution, each particle is supposed to have a diameter 
of 5 voxels. Due to reconstruction and scattering effects, 
the particles seem to be a little larger, with decreasing 
brightness towards the rim. However we can assume that 
the innermost voxel represents the gray value of the bulk 
material. This method leads to a volume filling factor 



of. 



0.37 



+0.06 
-0.05' 



in agreement with the method used in 



Paper I. While the volume filling factor is directly mea- 
surable, the coordination number was identified by Wada 
et al.| (I2OIII) to be the critical value for bouncing. Ac- 
cording to Eq. [l]our agglomerates correspond to a coor- 
dination number of Uc = 5.9, which is close to the critical 
value of Uc = 6. However, using the relations between 



spends to coordination numbers of ric = 3.76, ric = 4.92 
and Tic = 3.56, respectively. These values are signif- 
icantly lowe r than the critical value derived by |Wada 



et al. 



(2011). 



coordination number an d volume filling fact or by | Van de 



[Lagemaat et al.^ ( ,200 1] ), Rumpf ( p^958J and Yang et al. 



2.3. Inner Structure of the Dust Aggregates 

In Section [l] we discussed that the bouncing of dust 
agglomerates is currently debated between theorists and 
experimentalists. The question has risen whether bounc- 
ing, which has been observed in experiments but not in 
numerical collisions, is an artifact due to a compacted 
shell of the dust agglomerates. This compact mantle 
could be caused by the handling process or by the pro- 
duction method o f the dust aggregates as argured by 
Wada et al. (2011 ). We therefore analyzed the rim struc- 
ture of the nano-CT scanned dust aggregate at three dif- 
ferent positions. Figure [4] shows the density curve close 
to the surface and down to a depth of 18 /im, which 
corresponds to 12 monomer diameters. The density at 
the dust agglomerate's surface is close to the average 
value measured in the center of the dust agglomerate. 
Our measurements show a slight oscillation in the density 
profile. The variation in the volume filling factor is on 
the order of Acj) = ±0.05. With Eq. [ijthis corresponds 
to a deviation from the mean coordination number of 
Aric = ±0.4. We assume that the measured oscillation 
indicates a spatially regular structure close to the sur- 
face. 

2.4. Determination of the Collision Properties 

The main scope of our experiments was to study the 
outcome of collisions between sub-mm-sized dust ag- 
glomerates. We used a semi-automatic tracking pro- 
gram for individual dust agglomerates, which enables 
us to measure the collision velocities and dust-aggregate 
masses for identified collisions. However, only a small 
fraction of the actual number of collisions could be ob- 
served, due to optical-depth and identification problems. 
We are not aware of any preferred selection of a cer- 
tain collision type so that the observed collisions should 
be treated as a random representative (but incomplete) 
sample. 

After a collision was identified, the trajectories of the 
colliding dust aggregates were tracked before the colli- 
sion for as long as possible. From these trajectories, we 
calculated the collision velocities. It has to be taken into 
account that the dust agglomerates were only followed 
in two dimensions, due to the large number of individual 
dust agglomerates in both projected fields of view and the 
consequential impossibility to unambiguously identify a 
given d ust aggregate on the other projection. However, 
Weidling et al. (2012) showed that this restriction leads 
to a statistical error of only 13% in the collision velocity. 
We also neglected the rotation of the dust aggregates, 
which was clearly present but not measurable. We also 
checked the particle trajectories for any sign of accelera- 
tion or deceleration, due to gas drag or electric charges. 
All measured trajectories were linear in space and time 
within the accuracy of our measurements. 

The mass of the dust agglomerates was estimated from 
the average projected area when the particle was tracked. 
Assuming a spherical shape of the dust aggregates, this 
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Fig. 4. — Top: Measurement of the volume filling factor close to 
and below the surface of the analyzed dust agglomerate. The three 
curves represent the measurements for three different positions on 
the surface ( solid curve: 245 px x 265 px, dotted curve: 400 px x 
257 px, dash-dotted curve: 103 px x 142 px). The inset shows a 3D 
reconstruction of one of the analyzed regions. Bottom: NanoCT 
image of a cross-section through the dust agglomerate with semi- 
major axes of 150 /i,m, 100 /im, and 40 /im. 

leads to a dust- aggregate mass of 



M 



?SiO: 



^^3/2 



(2) 



where (j) = 0.37, ^si02 7 A are the volume filling factor 
of the dust aggregate, the mass density of the monomer 
grains, and the average projected area of the dust aggre- 
gate while the particle is tracked, respectively. 

3. RESULTS 

Figure |5] shows the temporal evolution of the four 
dust-aggregate dispersions during their respective experi- 
ments. The first column shows a snapshot after 1 second. 
The agglomerates made of polydisperse dust (Expl and 
Exp2) are well distributed throughout the experiment 
volume. In contrast, the agglomerates made of monodis- 
perse dust (Exp3 and Exp4) are only partially deagglom- 
erated. A large fraction of the mass in Exp3 and Exp4 
is still bound within a few large clusters, which we were 
not able to deagglomerate with our setup and the given 
experiment time. We call an agglomerate a cluster if it 
consists out of several agglomerates (like the large ag- 
glomerates in Fig. IgI) from the initially sieved particles 
(see Fig 11 d). 



As we will show in Section 3.3, individual particles 
are more likely to stick to clusters than to one another. 
Therefore, the single particles in Exp3 and Exp4 get de- 
pleted within a short time after the shaking frequency 
and, thus, also the average velocity has been reduced. In 
Exp3 we were able to observe the formation of a large and 
elongated cluster out of smaller dust agglomerates and 
clusters. This object has fractal-like characteristics and 
was formed in a succession of collisions as shown in Fig. 
[6] On the other hand, the polydisperse dust particles did 
not form very large clusters within the experiment du- 
ration. A reason for this might be that the shaking fre- 
quency was not reduced below 50% of the initial value. 
This caused the optical depth in these experiments to 
be higher than in Exp3 and Exp4 and made it difficult 
to identify individual collisions between dust agglomer- 
ates consisting of irregular Si02, albeit there should be 
many. For the experiments with monodisperse dust, the 
number of detectable collisions between individual sub- 
mm-sized dust aggregates is also limited, due to their 
low number density and due to the presence of the large 
clusters, which collide with the remaining particles more 
frequently du e to their larger cross section. 

In Section |3.1[ we will describe the outcome of the 
observed collisions between pairs of individual dust ag- 
glomerates. The properties of the cluster structure from 
a sequence of sticking collisions will be studied in Sec- 
tion |3.2[ A comparison between collisions of individual 
dust aggregates and tho se o f clusters of aggregates will 
be performed in Section [33] 



3.1. Sticking and Bouncing of Single Sub-mm-sized 
Dust Aggregates 

We analyzed a total number of 42 collisions between 
individual sub-mm-sized dust aggregates. The collision 
velocities were between 2.3 cm s~^ and 16.6 cm s~^ and 
the dust- aggregate masses were between 3.0 • 10~^g and 
9.3 • 10~^g. The collision partner had approximately the 
same size, with more than 85% of the particles having a 
mass ratio smaller than 5. The outcome of these colli- 
sions is presented in Fig. [7| in which the green and yellow 
data points denote sticking and bouncing collisions, re- 
spectively, and the long-dashed and short-dashed lines 
represent the 50%/50% and 0%/100% sticking/bouncing 
probabilities from Paper I, of which the parameters are 
given in Table |2j Triangular and circular symbols repre- 
sent collisions among dust aggregates consisting of poly- 
disperse and monodisperse Si02, respectively. A collision 
was identified as sticking when the particles sustained 
their contact for several frames. In particular, individual 
particle rotation after the collision illustrates that the 
sticking was impermanent, whereas a rotation of both 
colliding aggregates about a joint axis for an extended 
period of time indicates sticking. The observation time 
after a collision was limited by either a collision with a 
wall, a cluster covering the sticking agglomerates or the 
end of the experiment. The upper cloud of data points in 
the mass range around 10~^ g was taken from Paper I. As 
described earlier, our experiments were performed with 
two different dust samples. However, we found no strong 
evidence that the transition between sticking and bounc- 
ing depends on the monomer size distribution and mor- 
phology which is contrary to our observation that larger 
agglomerates only form in experiments conducted with 
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Fig. 5. — Snapshots of the four experiments at different times. The first two columns show the particles during the deagglomeration 
phase (100% shaker speed). The speed of the shaker is indicated by bars next to the images. For Expl and Exp2, the shaking frequency 
was reduced to 50% of the initial value after 5 seconds. In Exp3 and Exp4, the frequency was reduced to 50% of the initial value after 2 
seconds, to 35% after 5 seconds, and set to after 8 seconds. The right column shows the experiments at their final states. 
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Fig. 7. — Outcomes of the collisions between individual sub-mm- 
sized dust aggregates. The axes denote to the collision velocity and 
the mass of the smaller collision partner. Triangular and circular 
symbols represent collisions between dust aggregates consisting of 
polydisperse and monodisperse Si02 grains, respectively. The up- 
per cloud of data points are the collisions studied in Paper I. The 
long-dashed and short-dashed lines represent the 50%/50% and 
0%/100% sticking/bouncing probabilities from Paper I (see Table 
|2|. A new power law at 50%/50% sticking/bouncing probability 
was calculated to divide collisions which preferably lead to sticking 
for lower and bouncing for higher velocities. The fit is based on the 
results from this work (low masses) and Paper I (high masses). The 
dotted and solid lines indicate the continuous transition between 
perfect sticking (100%) and perfect bouncing (0%). 

the monodisperse dust partricles. This wih be diskussed 
in Section [4j It is obvious from Fig. [7| that cohisions 
among sub-mm-sized dust aggregates lead to sticking at 
velocities higher than predicted by Eq. 8 in Paper 1. 
Sticking also occurred at velocities for which pure bounc- 
ing was expected. 

3.2. Structures of Large Dust- Agglomerate Clusters 

In Section |3.1[ we showed that many of the observed 
collisions lead to sticking in the velocity range investi- 
gated. However, the final result of the growth process 
at the end of the experiment can be different. In Exp3 
and Exp4, we found that at the end of the experiment 
almost all initial dust aggregates were bound in large 
clusters (some of these clusters had survived the deag- 
glomeration phase at the beginning of the experiments). 
We also observed the formation of elongated and highly 
porous structures (see Fig. |6|. The formation of these 
fractal-like dust agglomerates has also been observed in 
later drop-tower experiments (Brisset et al., in prep). 
Measuring the fractal dimension of these large clusters is 
difficult, due to the finite initial dust-aggregate size and 
the limited resolution of the images. A rough estima- 
tion of the value of the fractal dimension can be derived 
by the box-counting method ( [Falconer 1990). For this 
method, the projected image ot the particle is placed in- 
side a two-dimensional box of the length 1. This box is 
divided in e grid cells with length 1/e. The number of 
grid cells is increased step by step. In each calculation 
step, the number N of boxes which cover a part of the ag- 
glomerate is counted. The fractal dimension is calculated 
from the relation N(e) = e~^f or taken from the slope of 
the linear function in a log (N) / log (1/e) plot. However, 
this method should only be treated as a rough estima- 
tion. The finite size of the ~ 150 /im agglomerates limits 



the resolution of the box-counting method. For Df < 2 
the value obtained from a projection of an agglomerate 
is comparable (but in some case s smaller) to the real (i.e. 
three-dimensional) value ( Nelson et aL J 1 9j0" ) . In our case 
the finite size of the initial dust aggregates requires an 
additional estimation of the error of our method. We 
applied the box-counting-method to every image of the 
agglomerate, which corresponds to images from different 
angles due to the rotation of the aggregates. Therefore, 
the value of the fractal dimension changes from image to 
image. These different perspectives allow us an estima- 
tion of the error of the method. The average value in 
the right image is Df = 1.58 ± 0.03 (more "side views") 
and Df = 1.71 ± 0.02 in the left image (more "head-on 
views"). The average of all images from both frames, 
considering the different number of available images in 
both projections, leads to a value Df = 1.63 ± 0.07. It 
must be noticed that there is also a systematic error due 



to the two-dimensional projection. Nelson et al. (1990) 



calculated that the deviation depends on the observation 
scale (the ratio of the size of the building blocks and the 
size of t he fractal) and the real value of Df (Eq.l3 and 
Fig. 2 in Nelson et al. (1990)). The error is negligible for 
values close to 2 and can be on the order of ADf =0.1 for 
fractal dimensions of Df 0.1. If we assume that the size 
of a typical individual monomer-agglomerate is on the 
order of 5 px and the analyzed fractal of 100 px, a mea- 
sured fractal dimension of Df = 1.63 would correspond 
to a real value of Df = 1.73. For the more compacted 
clusters, this method loses accuracy. Irregularities with 
sizes on the order of an individual dust aggregates on the 
surface of the projected cluster lead to lower values than 
expected from the images. Average values vary between 
Df = 1.63 and Df = 1.82 with an error that can be on 
the order of ADf = 0.08, at which we would expect val- 
ues closer to 2. 

The low fractal dimension for the elongated structures is 
in agreement with hit-and-stick cluster- c l uster collisions 
fBlum et al.|| 
minik 2006! 
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1 very low value is an indication that the 
threshold velocity between sticking and bouncing of clus- 
ters of dust aggregates is (well) above the experimental 
velocities. 

3.3. Sticking and Bouncing of Clusters of Dust 
Aggregates 



In Section 3.1, we focused on collisions between indi- 



vidual sub-mm-sized dust agglomerates. However, due 
to the succession of sticking collisions, we expect the 
dust agglomerates in protoplanetary disks to possess a 
more complex morphology than spherical. Therefore, 
collisions between clusters of dust agglomerates, which 
formed during the runtime of the experiment, are also 
of interest. We measured the relative velocities for col- 
lisions between dust- agglomerate clusters as well as for 
impacts of single dust aggregates into larger clusters. 

In the case of clusters of dust agglomerates the mass 
determination is difficult. We estimated the mass of the 
fractal-like objects by summing up the masses of all dust 
aggregates that were captured during the growth pro- 
cess (see Fig. |6|. For each of the captured elongated 
clusters, we mooeled the structure from the different im- 
ages available, by approximating the shape of the cluster 
with ellipsoids. The size of these "building blocks" was 
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Fig. 8. — Outcomes in collisions between clusters of sub-mm- 
sized dust aggregates consisting of monodisperse Si02 grains. The 
masses are those of the smaller cluster or, in the case of the five 
data points below 10~^ g, individual sub-mm-sized dust aggre- 
gates. The dotted and solid lines are the same as in Fig. mand 
indicate the continuous transition between perfect sticking (1D0%) 
and perfect bouncing (0%) in collisions between individual sub- 
mm-sized dust aggregates. 

estimated from the available images of the cluster. The 
mass estimation is the sum of these ellipsoids. In the case 
of the more compact clusters, we based the mass estima- 
tion on the projected area of the clusters and calculated 
the mass by Eq. [2] We likely overestimated the masses 
of these objects, but should still be on the right order of 
magnitude. 

The collision outcomes of the cluster impacts of Exp3 
and Exp4 are shown in Fig. [8] in the same way as in 
Fig. [7[ The green and yellow data points indicate stick- 
ing and bouncing collisions, respectively. The ordinate 
denotes the mass of the smaller agglomerate or cluster. 
The mass range of the larger collision partners was be- 
tween 4.4 • 10~^g and 1.3 • 10~^g. The five collisions with 
dust-agglomerate masses below 10~^ g are impacts of in- 
dividual dust agglomerates into clusters of agglomerates. 
All other impacts are between clusters of dust agglomer- 
ates. 

Compared with the data of the collisions among single 
dust aggregates (see Fig. |7|), we find many more sticking 
events than predicted by trie 50%/50% sticking/bouncing 
line (solid line in Fig. Isl). Within the range of our exper- 
imental data we were not able to see a clear onset of a 
sticking / bouncing transition for clusters. It is obvious 
that the looser structure of the clusters favors sticking, 
due to a better energy dissipation. 

4. DISCUSSION 

The experiments in this paper were performed with 
dust aggregates smaller, but structurally comparable to 
those used in Paper I. We were able to confirm a transi- 
tion between sticking and bouncing over a velocity range. 
However, a comparison with the transition curve and 
width in Paper I showed that their values (long and 
short-dashed lines in Fig. [t]) do not agree with the re- 
sults of this work. This disagreement is caused by the 
fact that the dust aggregates in Paper I had a rather 
narrow size distribution so that a deduction of the mass- 
velocity relation was not feasible. Thus, Pa per I used the 
contact model by jThornton fc Ning (1998) and its slope. 



which obviously is in contradiction with our new data. 

Now, we are in the situation that the combination of 
the data of Paper I and the data presented here spans 
more than three orders of magnitude in mass so that we 
can try and fit a power-law relation between the dust- 
aggregate mass m and the transition velocity v from 
sticking to bouncing, i.e.. 



m 



= a ■ 



1 cm/s 



(3) 



with a and b being free fit parameters. Due to the qual- 
itatively different evolution of the experiments with dif- 
ferent morphology of the individual dust grains (see Fig. 
[5|, we additionally consider a subset of data from ex- 
periments conducted with irregular dust only. We per- 
formed the fitting procedure with the logarithms of mass 
and velocity and considered the transition velocity as the 
50%/50% sticking/bouncing threshold. To calculate the 
power law, we used six different methods. For all fits, 
we minimized the quadratic and linear deviation as well 
as the number of data points on the wrong side of the 
fit function. For those three approaches, we minimized 
the error with and without the boundary condition that 
the number of false data points on both sides of the fit 
function is equal. 

In Table [2] we give an overview of the results of the dif- 
ferent methods. Values in parentheses were calculated by 
using only data from experiments with irregular dust par- 
ticles. To interpret these results, we used the conditional 
value at risk (CVaR) or expected shortfall (Hull 2007), 
which is often used in finance mathematics as a mea- 
sure to evaluate the risk of a portfolio. In our case, the 
CVaR is defined as the mean deviation in velocity of the 
q% collisions which deviate the most from the 50% line. 
Due to the small number of collisions deviating from our 
power laws, we used q=20%. Based on this criterion, the 
least squares and least linear deviation method and the 
requirement of a symmetric number of false data points 
on both sides of the curve as well as the least squares 
method with asymmetric false data points possess the 
lowest CVaR. The least-square and least-linear methods 
give rather similar exponents {b ~ —4/3) and prefactors 
(a ~ —4.16... — 4.79) which are, within the accuracy of 
our measurements and taking the small number of data 
points into account, in good agreement. Some of the cal- 
culation methods lead to several solutions within small 
derivations in a and b. For these cases we give an aver- 
age value of the parameters and marked them with an 
asterisk in Tabled 

We also calculated the threshold power law for the col- 
lisions from Paper I and only the results from the exper- 
iments with polydisperse dust. This leads to comparable 
values for the prefactors a. However, the power-laws dif- 
fer slightly from the values obtained with the full data 
set but do not change the overall picture (see Table 

In the same fashion as in Paper I, we assume a log- 
arithmic sticking-probability distribution in velocity so 
that all threshold lines are parallel. However, we were 
not able to treat the dust-aggregate masses ranging over 
more than three orders of magnitude as identical. For 
each sticking and bouncing collision, we therefore calcu- 
lated the deviation in velocity according to the power law 
computed above. As shown in Fig. [9) we then used log- 
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TABLE 2 

Fit results of the six different methods used to calculate the transition between bouncing and sticking for monodisperse 
and polydispers dust samples. values in parentheses are only based on the polydisperse sample from this work and 
[Weidling et al.] ( |2Q12[ ). Calculations marked with asterisks represent the average of several similar solutions. For 

details, refer to the text. 



Fit 


Exponent 


Constant 


Tran- 


CVaR 


method 


b 


a 


sition 
width r 


[log(cms-i)] 


Asymmetric number of false data points 


Least squares 
Least linear dev. 


-1.34 (-1.13) 
-1.36 (-1.13) 


10-4.34 (10-4.79) 
10-4.42 (10-4.79 ) 


-0.72 (-0.44) 
-0.78 (-0.45) 


0.41 (0.86) 
0.46 (0.88) 


Least number of 


-0.94 (-0.93) 


10-4-87 (10-4-88 ) 


-0.40 (-0.38) 


1.12 (1.15) 


false data points* 










Symmetric number of false data points 


Least squares* 
Least linear dev.* 
least number of 


-1.35 (-1.48) 
-1.37 (-1.47) 
-1.30 (0.91) 


10-4-31 (10-4-33) 
10-4-31 (10-4-16 ) 
10-4.46 (10-4.30 ) 


-0.71 (-0.71) 
-0.73 (-0.69) 
-0.80 (-0.43) 


0.39 (0.41) 
0.39 (0.45) 
0.51 (0.72) 


false data points* 










Values from Paper I 




-3.6 


10-7.76 







arithmic equidistant velocity bins to fit a linear function 
for the probability distribution. 



P = r • A + 0.5, 



(4) 



where P is the sticking probability between and 1, A 
is the logarithmic deviation in velocity from the P=0.5 
power-law and r = —0.73 the slope of the function (least 
linear deviation and symmetric fit method). The tran- 
sition function P(A) is plotted in Fig. ItI The slight 
disagreement between the data points anathe fit curve 
is caused by the fixed value for the 50%/50% stick- 
ing/bouncing threshold. For comparison we also show 
the sticking probability, within the velocity bins, for 
our experiments with only monodisperse (green dots) 
and polydisperse (red diamonds) dust. Considering the 
smaller number of collisions, these data subsets are still 
described by the fit derived for all data, even though 
the agglomerates made out of monodisperse dust stick 
at slightly higher velocities. 

The new power law for the 50%/50% sticking/bouncing 
transition is less steep (slope ^ ~4/ 3 ) than the one given 
by the model of Thornton fc Ning] ([19981), with a slope 
of -18/5, used in Paper 1. Small dust aggregates stick at 
higher velocities than predicted by Paper I. This shows 
the limits of the previous contact model, which is based 
on static contacts between solid spheres. Moreover, the 
transition from perfect sticking to bouncing is much nar- 
rower and the reason is the lack of data for sticking col- 
lisions in Paper I. It is obvious for the two clouds of data 
that the transition can be better determined when in- 
cluding the small- aggregate collisions. 

The impact of this new transition curve between stick- 
ing and bouncing and the much narrower transition 
width on our understanding of the formation of small 
bodies in the solar nebular is difficult to predict. Also two 
potential limitations of our results should be considered. 
The experiments were conducted with two different dust 
analog materials, both consisting of silica but one with 
a monodisperse and one with a wider size distribution. 
While the exact material is supposed to have a negligible 
impact on the experiment outcome (e.g. Blum & Wurm 
f2QQ8J ), the size distribution and the monomer size and 
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Fig. 9. — Calculation of the sticking probability based on the 
data from this work and Paper I. The abscissa shows the logarithm 
of the deviation from the best fitting power law (least linear dev., 
symmetric) for the 50%/50% sticking/bouncing threshold. Plot- 
ted are the results of the individual collisions (black bars) and the 
sticking probability within equidistant logarithmic velocity bins for 
both materials (blue) and only the monodisperse (green dots) and 
polydisperse (red diamonds) from this work. To describe the tran- 
sition from sticking to bouncing, we applied the linear function 
shown in Eq. [4] The 50%/50% sticking/bouncing position was 
fixed by Eq. [3] 

shape clearly have. It is assumed that the dust pop- 
ulation in the protoplanetary disks consists of particles 
similar to today's interstellar dust grains, with sizes f rom 
the nanome ter range up to about 1 micrometer. (Ma this| 



et al. 



1977). However, the abundance of even larger 



grams is not well determined. Studies o f impacts of in- 
terste llar dust grains in the solar system fLand^raf et al 



2000 ), of primitive, chondritic meteorites (Scott &; Krot 



2005) and the samples returned from comet 81P / Wildz 
by the Stardust mission (Brownlee et al. 2006) showed 
larger dust grains, up to a size ot ^ 10/im. Assuming 
that these particles are the remains of the initial proto- 
planetary dust population, we can assume grain sizes not 
entirely different from the ones used in our experiment. 
An estimation of the influence of the monomer size on 
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the sticking velocity is difficult. In principle, a larger 
diameter corresponding to the average monomer mass 
of the polydisperse analog dust (2.0 fim) should lead 
to a lower maximum sticking velocity compared to the 
monodisperse dust (1.5 /am) (compare e.g. the sticking / 
fragmentation transition in [Dominik fc Tielens (1997)). 
However, due to the limited data it is not possible to see 
a clear trend in our data. Figure [9] shows that the ag- 
glomerates made out of polydisperse dust (red diamonds) 
stick at slightly lower velocities than the agglomerates 
made of monodisperse particles (green dots), which is in 
agreement with theory. However, a model-based-, quan- 
titative comparison of the two threshold velocities is, un- 
fortunately, not feasible, as theoretical models assume 
spherical dust grains. This is not the case for our poly- 
disperse dust samples which consist of irregular shaped 
grains. 

A second point which should be addressed is the single 
porosity of our samples ^ 0.37) . Various numeri- 
cal simulations ([Suyama etal. 2008 Zsom et al. 2011 



Okuzumi et al. 



2012p have shown that the first dust ag- 



glomerates might have had a much lower volume filling 
factor (down to (j) ~ 10~^). Zsom et al.| ( [201 1 ) showed 
that particles with a mass of ^ 10~^ g can still have 
such a low density. Their simulations also showed that, 
after a relatively short time (few hundred years), these 
particles can become more compact and end with val- 
ues comparable to our samples. However, this is based 
on the compac t ion m odel from Guttle r et ah] ( |2010t ) and 



Weidling et al.| (|2009 ) who already assume spherical com 
pacted [(j) ^O.i) particles. Our experiments showed that 
even aggregates of a few hundred micrometers in size are 
more likely to form open structures than spherical ob- 
jects and thus have a lower filling factor than previously 
assumed. The formation of very small, fractal agglomer- 
ates had also been obs erved in experiments studying an 
earlier growth regime ( iWurm fc Bluni|[T998l [Blum et al. 



2002). In our case, collisions of aggregate clusters, which 
led to sticking in most cases, caused hardly any visible 
compaction. This fact emphasizes that also small ag- 
glomerates might have had a very low filling factor for a 
longer time than predicted in simulations, because they 
were not in a regime in which bouncing would compact 
them. Therefore, it should be considered that our new 
collision model might only be applicable to later phases 
of planetesimal formation after the dust aggregates got 
compacted. 

Considering these constraints, our results suggest that 
there is a region of increased collisional growth between 
equal-sized dust aggregates which benefits the forma- 
tion of larger objects. When dust aggregates possess 
a collision-velocity distribution rather than a fixed con- 
stant value, a combination with the probability distri- 
bution of sticking and bouncing can lead to a contin- 
uous growth of " l ucky w inners", as recently shown by 
Windmark et al.| (|2012b|). Dust-aggregate growth can 



also benefit from the open structure of the growing dust 
agglomerates. Current models make the assumption that 
a dust aggregate of a given mass has a spherical shape 
and homogenous structure. An open-structured cluster 
of dust aggregates of the same mass would more likely 
stick in a collision due to its ability to dissi pate more 
energy in restructuring, as shown in Section |3.3| Such 



clusters also possess a larger cross section, which makes 
collisions more likely. This effect might be of interest for 
future simulations. However, more studies of the fractal 
dimension and the restructuring of these clusters of dust 
agglomerates are needed. Fig. [Tjalso indicates that there 
might be a direct transition from sticking to fragmenta- 
tion for dust aggregates with masses below ~ 10~^ g. 
This mass corresponds to a dust agglomerate which con- 
sists of ^ 10^ monomer grains. Such a dust agglomerate 
is about ten times as massive as the ones simulated by 
Wada et all (|2011|), which can possibly explain why they 
did not observe bouncing for porous dust agglomerates. 
Anyhow, it is questionable if the onset of fragmenta tion is 
totally mass independent. Recen t experiments by [Beitz 
et al.| (|2011') and Schrapler et al.| ([2Q12) have shown that 
dust aggregates larger than one centimeter fragment at 
velocities slightly lower than 1 ms~^. This, however, 
does not change the picture qualitatively. 

5. CONCLUSION 

We measured the outcome of collisions between sub- 
mm-sized dust agglomerates at low velocities of sev- 
eral cm s~^ under micro-gravity conditions. We were 
able to confirm the findings of earlier experiments with 
larger dust aggregates that a transition region between 
perfect sticking and perfect bouncing exists. Based on 
our experiments and those from Paper I, we determined 
a new power law for the mass-velocity dependence of 
the 50%/50% sticking/bouncing value as well as the 
width of the transition from 100% /0% to 0%/100% stick- 
ing/bouncing. This power law is shallower than assumed 
in Pap er I which is based on the model by |Giittler et al. 
( |2010| . This suggests that growth due to sticking is pos- 
sible even at velocities higher than before assumed. This 
disagreement also shows that the collision physics be- 
tween porous agglomerates can not be easily described by 
modified contact models for solid spheres. We also found 
that clusters of dust aggregates even stick at each other 
at velocities where compact dust agglomerates would 
only bounce. Previous collision models, which were used 
in growth simulations, treated the newly formed agglom- 
erates as compact spheres. Considering our observa- 
tion that open structures benefit the growth significantly 
might help to overcome the various barriers found in sim- 
ulations. These clusters can either be rather spherical or 
open-structured. For the latter species, we analyzed the 
morphology of one example and found that its fractal 
dimension is on the order of Df = 1.7. Furthermore, we 
analyzed the structure of one of our dust aggregates using 
nano-CT imaging. We confirmed previous measurements 
of the filling factor to be about (j) = 0.37 and were not 
able to find any compacted rim that had been suggested 
to be a reason for the bouncing of dust aggregates. We 
also estimated the coordination number of these dust ag- 
gregates. Using different relations from the literature, it 
is most likely that our dust aggregates possess a coordi- 
nation number well below Uc = 6, which is assumed to be 
the critical value at which simulations showed bouncing. 
Our new data for the sticking-bouncing threshold suggest 
that a direct transition from sticking to fragmentation, 
as found in numerical simulations, is possible for dust ag- 
gregates consisting of less than ~ 10^ monomer grains, 
which is the range of the numerical simulations. Thus, 
the non-occurrence of bouncing could be a finite-size ef- 



13 



feet. 
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